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ABSTRACT 

We use the supernova remnants (SNRs) in the two Magellanic Clouds (MCs) as a 
supernova (SN) survey, "conducted" over tens of kyr, from which we derive the current 
SN rate, and the SN delay time distribution (DTD), i.e., the SN rate vs. time that 
would follow a hypothetical brief burst of a star formation. In a companion paper 
(Badenes, Maoz, & Draine 2010) we have compiled a list of 77 SNRs in the MCs, and 
argued that it is a fairly complete record of the SNRs that are now in the Sedov phase 
of their expansions. We recover the SN DTD by comparing the numbers of SNRs 
observed in small individual "cells" in these galaxies to the star-formation histories 
of each cell, as calculated from resolved stellar populations by Harris & Zaritsky. We 
identify the visibility times of SNRs in each cell with the Sedov-phase lifetimes, which 
depend on the local ambient densities. The local densities are estimated from 21-cm 
emission, from an inverse Schmidt-Kennicutt law based on either Ha emission or on 
the star-formation rate from the resolved stellar populations, and from combinations of 
these tracers. This is the first SN DTD that is based on resolved stellar populations. 
We detect a population of "prompt" type-la SNe (that explode within 330 Myr of 
star formation) at > 99% confidence level (c.L). The best fit for the number of prompt 
type-la SNe per stellar mass formed is 2.7— 11.0 x 10 _3 Mq , depending on the density 
tracer used. The 95% c.l. range for a "delayed" (from 330 Myr to a Hubble time) type- 
la component is < 1.6 x 10~ 13 SN yr _1 MQ , consistent with rate measurements in old 
populations. The current total (core-collapse+Ia) SN rate in the MCs is 2.5-4.6 SNe 
per millenium (68% c.l.+systematics), or 1.7-3.1 SNuM [SNe (100 yr IO^Mq)" 1 ], in 
agreement with the historical record and with rates measured in other dwarf irregulars. 
Conversely, assuming the SNRs are in free-expanion, rather than in their Sedov phase, 
would impose on the SNRs a maximum age of 6 kyr, and would imply a MC SN rate 
per unit mass that is 5 times higher than in any type of galaxy, and a low-mass limit 
for core-collapse progenitors in conflict with stellar evolution theory. 

Key words: supernovae: general - supernovae: remnants galaxics:individual: LMC, 
SMC 



1 INTRODUCTION 

Supernova (SN) explosions and their remnants touch upon 
multiple aspects of astrophysics and cosmology, whether as 
endpoints in stellar evolution, as sources of energy and en- 
riched material to the interstellar and intergalactic media, 
as sites of cosmic-ray acceleration, or as standard candles 
for cosmography. However, much remains to be understood 
regarding these events, both the core-collapse (CC) SN ex- 
plosions that are thought to end the lives of some or all 
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massive (> 8Mq) stars; and the type-la SNe (SNe la), 
that are believed to be the thermonuclear combustions of 
CO white dwarfs (WDs) that have approached the Chan- 
drasekhar mass by accreting material from, or merging with, 
a companion star. 

The lower limit on the initial mass of stars that even- 
tually undergo CC is poorly known, both theoretically 
l|Poelarends et ail 120081 ) and observationally (from identi- 
ficati on of progenitors in pre-explosion images; e.g. ISmartj 
I2009D . and could be as low as 7 Mq or as high as 12 Mq . The 
upper mass limit on the intial mass that leads to a SN ex- 
plosion, rather than (perhaps) to an explosion-less collapse 
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into a black hole, is also unc ertain theoretically |fleger et alj 
I2003T ) and observationally (|Kochanek et alj|2008h . Among 
the CC-SN progenitors, it is still not clear which mass range 
leads to whic h SN subtypes, such as IIP, Iln, lb, and Ic 
(|Smarttll2009h . Parameters other than mass - binarity, ro- 
tation, metallicity - likely al so play an importan t role in 
determining CC-SN type (e.g. lEldridge et alJliooih . 

For SNe la, our ignorance is even greater. The pro- 
genitor systems of these events are still unknown, with two 
distinct models generally considered: the single-degenerate 
(SD) scenari o in whi ch the WD accretes material from 
normal star l|Whelan 



Table 1. Stellar Lifetimes 



Ibenlll973l ; iNomotol Il982r i 



a 

and the 

double-degene rate (DP) model, where the WD merg es with 
another WD (|lben fc Tutukovlll984l ; IWebbinklll984r ). Both 
models suffer from theoretical and observational problems. 
For SD scenarios, it is unclear how the material accreted 
from the donor can burn quietly on the WD surface un- 
til the req uired mass (close t o the Chandrasekhar limit) 
is reached (|Cassisi et al.lfl99"sl ). The narrow H or He lines 
that are expected in nebular SN spectra for SD progeni- 
tors have never been found in normal la events (I Leonard! 
I2007D , nor have the integrated X-ray emission from a pop- 
ulation of accreting WDs undergoing slow nuclear burn- 



ing o n their surfaces (|Di Stefanoll2010l ; iGilfanov fc Bogdan 
120101 ). The claimed identification of the surviving donor 
star in the nearby Type la supernova remnant (SNR) Ty- 



cho (|Ruiz-Lapuente et al . 2004; iHernandez et alIl2009T) has 



been recently questioned bv lKerzendorf et al. I i|2009h . In the 
case of DD systems, there are concerns regarding the fate of 
merging WDs, which might lead to accretion induced co l- 
lapse instead of a SN la explosion (jSaio fc Nomoto|[l985l ). 
It is also unknown whether there are enough DD progeni- 
tor candi dates in the Milky Way to produce the observe d la 
SN rate l|Napiwotzki et all [2004 ; iNelemans et alj|2005l ). al- 
though currerrt_se^xchesjrright clarify this point in the near 
future (|Badenes et al.ll2009bh . 

Major progress in resolving these questions could come 
from knowledge of the elapsed times between the formation 
of a stellar population and the explosion of some of its mem- 
bers as different types of SNe. Indeed, a major objective of 
SN studies has been the recovery of the so-called delay-time 
distribution (DTD). The DTD is the SN rate as a function 
of time that would be observed following a 5-function burst 
of star formation. (In other contexts, the DTD would be 
called the delay function, the transfer function, or Green's 
function). Knowledge of the DTD would be useful for un- 
derstanding the route along which cosmic metal enrichment 
and energy input by SNe proceed, but no less important, for 
obtaining clues about the SN progenitor systems. Different 
progenitor stars, binary systems, and binary evolution sce- 
narios predict different DTDs. The "border" in progenitor 
masses between SNe la and CC-SNe will lead to a time bor- 
der in the DTD, dictated by stellar evolution, between the 
two types. The age corresponding to a particular mass bor- 
der will depend on metallicity. For example, an 8M©, OAZq 
(where Zq is the solar metallicity) star lives for 41 Myr. 
For reference here and i n the rest of this paper, Table Q] lists 
the (|Girardi et al.ll2000h lifetimes of stars of various zero-age 
main sequence masses, for various metallicities. The mean 
metallicities of the Small and Large Magellanic Clouds, the 
two galaxies that will be at the focus of this paper, are O.IZq 
and OAZq, respectively (Russell & Dopita 1992). 
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The shape of the DTD is especially important to con- 
strain the hotly debated nature of SN la progenitors. For 
both of the currently popular progenitor scenarios, SD and 
DD, calculations of the DTD depend on a series of as- 
sumptions regarding initial conditions (initial mass func- 
tion - IMF , binarity fraction, mass ratio distribution, 
separation distribution), and complex physics (mass loss, 
mass transfer, common envelope evolution, accretion) that 
are often computationally intractable except in the most 
rudim e ntary, param et rized fo rms (e.g., Yungelson fc Livid 
19981; iHurlev et all [20021; lHan fc Podsiadlowskil I 206X 



Nelemans fc Toutll2005l ; iBeer et alj|2007l ; iRuiter et al.ll200S 
Bear fc Sokerl |2009| . Bogomazov & Tutukov 2009; Men- 
nekens et al. 2010). Observational estimates of the DTD 
could rule out particular DTD predictions, or at least could 
provide some input and generic features that successful mod- 
els will need to reproduce. 

However, to date, observational SN DTD estimates have 
been few and controversial. One approach has been to com- 
pare the SN rate in field galaxies, as a function of redshift, 
to the cosmic star formation history (SFH). Given that the 
DTD is the SN "response" to a short burst of star forma- 
tion, the SN rate versus cosmi c time will be the convol ution 
of the full SFH with the DTD. lGal-Yam fc Maozj (|2004h car- 
ried out the first such comparison, using a small sample of 
SNe la out to 2 = 0.8, and concluded that the results were 
strongly depende nt on the poor l y kno wn cosmic SFH, a con- 
clusion echoed bv lForster et al.l (|2006j ). Wit h the availability 
of SN rate measurements to higher redshift s . I B arris fc Tonrvl 
l|2004h found a SN la rate that closely tracks the SFH out to 
z ~ 1, and concluded that the DTD must be concentrated 
at short delays, < 1 Gyr. Similar conclusions have bee n 
reached, at least out to z ~ 0.7 , by Sullivan et al.l ( 2006). 
In con trast, iDahlen et all (|2004 l2008h and IStrolger et all 
(|2004h have argued for a DTD that is peaked at a delay of 
~ 3 Gyr, with little power at short dela ys, based on a de- 
crease in the SN la rate at z > 1. However. iKuznetsova et ahl 
(|2008h have re-analyzed some of these datasets and con- 
cluded that the small numbers of SNe and their potential 
classi fication errors p r eclud e reaching a conclusion. Simi- 
larly, IPoznanski et al.l (|2007f ) performed new measurements 
of the 2 > 1 SN la rate, and found that, within uncertainties, 
the SN rate could be trac king the SFH. This, again, would 
imply a short delay time. iGreggio et al.l (|2008T l pointed out 
that these results could also be affected by an underesti- 
mated extinction for the highest-2 SNe, which are observed 
in their rest-frame ultraviolet emission. 

A second approach to recovering the DTD has been 
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to compare the SN rates in stellar populations of differ- 



200dl200d ) 
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Scarmapicco & Bildstcn ( 20051) . ISullivan et al 



ent characte r istic ages. Using this approach. IMannucci et al 



2006), and lRaskin et al.l |2009|) have all found evidence for 



the co-existence of two SN la populations, a "prompt" pop- 
ulatiorQ that explodes within of order 10 8 yr, and a de- 
layed channel that produces SNe la on timescales of or- 
der 5 Gyr. Naturally, these two "channels" may in real- 
ity be just integrals ov er a continuous DTD on two sides 
of som e time border jGreggio et al.l 120081 ). iTotani et al.l 
|2008l ) have used a similar approach to recover the DTD, 
by comparing SN la rates in early-type galaxies of differ- 
ent characteristic ages. They find a DTD consistent with 
a t -1 form. Additional recent studies using this approach, 
but which may be influenced by selection effects and by 
a posteriori statistics (because they focus on the proper- 
t ie s of SN host g alaxie s ; see M aoz 2008) can be found 
in lAubourg et alj J200sl), ICooper et all (|2009l ). Schawinski 
(2009). and lYasuda fc Fukugital (|2010l ). 



A third approach for recovering the DTD is to measure 
the SN-rate vs. redshift in massive galaxy clusters. Opti- 
cal spectroscopy and muliwavelength photometry of cluster 
galaxies has shown that the bulk of their stars were formed 
withi n short episodes at z ~ 2 — 3 (e.g. lEisenhardt et all 
120081 ). Thus, the observed SN rate vs. cosmic time (since the 
stellar formation epoch) essentially provides a direct mea- 
surement of the form of the DTD. Furthermore, the record 
of metals stored in the intracluster medium (ICM) con- 
strains the number of SNe th at have exploded, and hen ce the 
normalization of the DTD. iMaoz fc Gal- Yarn! (|2004) anal- 



ysed the then-available cluster SN la rates ( Gal- Yam et al.l 
120021 ) . The low observed SN la rates out to 2 ~ 1 im- 
plied that the large number of events needed to produce 
the bulk of the iron occurred at even higher redshifts, be- 
yond the range of the then-existing observations. They con- 
cluded that most of the cluster SNe la exploded during the 
relatively brief time interval between star formation in mas- 
sive clusters (at z ~ 2 — 3) and the highest-redshift clus- 
ter SN rate measurements (at z ~ 1). Cluster SN rate 
measurements have by now greatly im proved in the red- 
shift range from zero to 1, and beyond JSharon et al.l 20071; 
Gal- Yam et all 120081; IMannucci et al.l l200Sl: iGraham et all 



20091 ; ISharonl l201(il ; iDildav et all 120101 ). Analysis of these 



rates (Maoz et al. 2010b) reinforces the conclusion that, al- 
though the SN la DTD may have a low-level tail at delays of 
a few Gyr, the bulk of the events must occur within ~ 1 Gyr 

of star form ation. 

Finally, IMaoz et al.1 (|2010al ) have presented a method 
for recovering the DTD from a SN survey, where the SFHs 
of the individual galaxies surveyed are taken into account. 
They applied the method to a subset of the Lick Obser- 
vatory SN Search (LOSS - Filippenko in prep.; Leaman 
et al. in prep.; Li et al. in pre p.) that overlap s with the 
Sloan Digital Sky Survey (SDSS; lYork et al.ll2000l ), for which 
SFH rec onstruction is avail ale based on the SDSS spec- 
troscopy (|Toieiro et alj|2009h . They found that a "prompt" 



(< 420 Myr) SN la component is required by these data at 
the > 99% confidence level. In addition, a delayed SN la 
population, with delays of > 2.5 Gyr, is detec ted at the 4a 
level. A related approach has been used by iBrandt et al.l 
(|2010l ). with similar conclusions. 

In this paper, we seek to recover the SN DTD in yet an- 
other way, by analyzing the SNRs in the Magellanic Clouds 
(MCs). The MC SNRs present several advantages. First, 
they are all at the known distances of their two host galax- 
ies. Second, both MCs have been surveyed to large depths 
in the radio, with individual sources followed up at multiple 
wavelengths and classified. In a companion paper (Badenes, 
Maoz, & Draine 2010, hereafter Paper I), we have compiled 
a sample of MC SNRs, and argued that it is largely com- 
plete. Third, the MCs are close enough to permit a region- 
by-region fitting of their resolved stellar populations with 
detailed stellar evolution isochrones. The SFHs of the Clouds 
can therefore be reconstructed with bett er spatial and tem- 
oral detai l than in any othe r galax ies. lHarris fc Zaritskvl 
2004! ) and lHarris fc Zaritskvl (|2009l ) have recently carried 
out such a program of SFH reconstruction of the MCs. 
iBadenes et"afl (|2009al ) have compared the SFHs of the MC 
regions hosting specific SNRs to the properties of the rem- 
nants, in order to deduce constraints on the nature of some 
of the explosions and their delay times. 

Here, we go a step further, and treat the MCs and 
their SNRs as an effective SN survey conducted in a sam- 
ple of galaxy subunits, where the detailed SFH of each 
subunit is known. In §2, bel ow, we briefly rev i ew th e re- 
construction of MC SFH s by lHarris fc Zaritskvl (|2004l ) and 
lHarris fc Zaritskvl l|2009l ). In §3 we review our MC SNR sam- 
ple from Paper I, and the physical model that reproduces 
the observed size distribution of this SNR sample. We show 
how this same model permits estimating the relative visibil- 
ity times of SNRs at different locations in the MCs, given the 
local ambient densities. As in Paper I, we use three different 
tracers of MC density to estimate the SNR visibility times 
as a function of lo cation. In §4, we review the method of 
IMaoz et all (|2010al ) for recovering the most likely SN DTD 
and its uncertainty, given the spatially resolved SFHs, the 
SNR visibility times, and the observed number of SNRs. In 
§5 we apply the method to the MCs with their sample of 
SNRs, and derive the DTD. In §6 we use the visibility times 
to obtain also the current SN rate in the MCs, and discuss 
the emerging picture. 

For the purpose of this paper, we assume that the DTD 
is a universal function: it is the same in all galaxies, inde- 
pendent of environment, metallicity, and cosmic time — a 
simplifying assumption that may be invalid at some level. 
For example, a dependence of SN delay time on metallicity 
is expected in some models (e.g., Kobayashi et al. 2000). 
Similarly, variations in the initial mass function (IMF) with 
cosmic times or environment would also lead to a variable 
DTD, but we will again ignore this possibility in the present 
context. 



We note that different authors have used the term "prompt" to 
describe different delays, from < 100 Myr to as long as < 1 Gyr. 
In this work, we will refer as prompt SNe la to those with delays 
< 330 Myr. 



2 THE STAR-FORMATION HISTORY MAPS 
OF THE MAGELLANIC CLOUDS 

The SFH maps that we use in the presen t wor k are 
presented in full detail in lHarris fc Zaritskvl l|2004l ) and 
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lHarris fc Zaritskvl :2(H)!lQ. The maps were elaborated us- 
ing four-band ( U, B, V, and /) phot ometry from the Mag - 
ellanic Clouds Photometric Survey l)Zaritskv et all 12004 ). 
which has a limiting magnitude between 20 and 21 in 
V, depending on the local degree of crowding in the im- 
ages. In each Cloud, the data were divided into regions or 
"cells" with enough stars to produce color-magnitude di- 
agrams of suffic ient quality, which wer e then fed into the 
StarFISH code (|Harris fc Zaritsky|[200ll ') to derive the local 
SFH for each cell. For th e Large Magellanic Cloud (LMC), 
lHarris fc Zaritskvl (|2009f ) divided more than 20 million stars 
into spatial cells encompassing the central 8° x 8° of the 
galaxy (see their figure 4). The majority of the cells are 
12' x 12' squares, while about 50 cells in regions of lower 
stellar density have sizes of 24' x 24'. In total, there are 1376 
cells for the LMC, for which the SFH is given in 13 tempo- 
ral bins with lookback times between 6.3 Myr and 15.8 Gyr. 
For th e Small Magellanic Cloud (SMO. lHarris fc Zaritskvl 
(|2004l ) divided over 6 million stars into 351 12' x 12' cells, 
leaving out two areas that are contaminated by Galactic 
globular clusters in the foreground (see their figure 3). The 
temporal binning of the SFHs is slightly different for the 
SMC, with 18 resolution elements between 4.6 Myr and 9.7 
Gyr. In the original maps, the SFHs are further subdivided 
into metallicity bins, but we ignore this dimension, working 
with the co-added SFHs instead. 

The stellar masses f ormed in every t i me b in in 
every cell, as de r ived by lHarris fc Zaritskvl (120041 ') and 

IMF. 



assume a 



Salpeterl 



lHarris fc Zaritskvl (|2009l ). 
For ease of comparison of our results to other SN rate work, 
we co nvert these masses to a "diet Salpeter" IMF (jBell et alj 
l2003h . by multiplying the masses formed by factor of 0.7, to 
account for the reduced number of low-mass stars i n a real- 
istic IMF, compared to the original ISalpeterl (| 19551 ) IMF. 



3 THE SAMPLE OF SUPERNOVA 

REMNANTS IN THE MAGELLANIC 
CLOUDS, ITS PROPERTIES, AND PHYSICS 

In Paper I, we assembled a list of 77 SNRs in the MCs, 
by cross- identifying literature compilations of SNRs. We 
showed that there is an observed "floor" of 50 mjy in the 
radio flux of the SNRs, orders of magnitude higher than the 
flux limits of the radio surveys in which almost all of these 
SNRs are detected. From this, we concluded that the paucity 
of SNRs below ~ 50 mjy must be real, rather than being 
due to any observational effect. This suggested that any MC 
SNRs with radio fluxes under the ~ 50 mjy floor but above 
the detection limit must fade very quickly through the flux 
range from the floor down to the detection limit. 

Analysing our sample in Paper I, we showed that the 
size distribution of the SNRs is close to linear in the cumula- 
tive, or uniform in the differential, up to a cutoff at a radius 
of r ~ 30 pc. Such a distribution of SNR sizes has been 
noted before both in the MCs and in other galaxies. It has 
been variously attributed to observational selection effects 
or to the SNRs being in a "free-expansion" phase (i.e., with 
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constant shock velocity), expected when the mass swept up 
by the expanding shock is still small compared to the ejected 
mass. We suggested an alternative physical model to explain 
the MC SNR size distribution. 

Briefly, we proposed that most of the MC SNRs are 
in the decelerating Sedov- Taylor phase, where the swept- 
up mass is larger than the ejecta mass, but the cooling 
time of the shocked gas is still longer than the age of the 
remnant, and hence the evolution is approximately adia- 
batic. Once the cooling time becomes comparable to the 
age, the SNR quickly loses energy radiatively, entering the 
radiation-loss-dominated snowplough phase, after which it 
slows down, breaks up, and merges with the interstellar 
medium. Quantitative estimates show that typical SNRs 
should be in their Sedov- Taylor phases for ages between a 
few and a few tens of kyrs , and for sizes of order a few to 
a few tens of parsecs fe e. ICioffi et al J 1 19881 ; iBlondin et al.l 
ll998tlTruelove fc McKeell 19991 ). During this phase, the SNR 
sizes grow as 



pl/5 -1/5,2/5 



(1) 



where Eo is the kinetic energy of the explosion, p is the 
ambient gas density, and t is the time. The shock velocity 
therefore decreases as 



dr 
~dt 



pl/5 -1/5,-3/5 



or equivalently expressed in terms of r rather than t, 

,-,1/2 -1/2 -3/2 

v ~ E ' p ' r ' . 



(2) 



(3) 



We proposed in Paper I that the uniform size distribu- 
tion arises as a result of the transition from the Sedov phase 
to the radiative phase. The radius of this transition depends 
on ambient density. When coupled with the distribution of 
densities in the MCs, this leads to fewer and fewer sites at 
which large-radius Sedov-phase SNRs can exist. The cooling 
time of the shocked gas depends on the density as 



tc 



kT 
pA(T) 



(4) 



where A(T) is the cooling function at temperature T. The 
cooling function can be approximated with a power law, T £ , 
in the temperature range of relevance for the shocked gas, 
around 10 6 K. By relating the temperature to the shock 
velocity v, kT ~ m p v 2 (where m p is the proton mass), and 
equating t coo \ to the age t as expressed in Eq. [2] one obtains 
that the transition radius, r max , scales as 



p (3-2e)/(ll-6 e ) -(5-2e)/(ll-6e) 



(5) 



at 



A remnant expanding beyond this radius, at the given am- 
bient density, will enter the radiative phase and quickly fade 
from view. For a fairly large range of plausible cooling func- 
tion dependences, e.g., indices e of —1/2 to —3/2, r max de- 
pends on density as p~ 3 ^ 7 to p~ 2 ^ 5 . Conversely, the maxi- 
mum ambient density that will permit a Sedov-phase SNR 
of radius r is 

Pmax oc r s , (6) 

where S is likely in the range —7/3 to —5/2. We showed in 
Paper I that a uniform size distribution, dN/dr ~ const., 
is obtained if the ambient gas density follows a power-law 
distribution of density, with index /3 ~ — 1. 
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We then examined three different tracers of gas density 
in the MCs: the HI column density; the star-formation rate 
(SFR) based on resolved stellar populations, translated to 
a density via an inverse Schmidt law; and the SFR based 
on the Ha emission-line luminosity, again combined with an 
inverse Schmidt law. We showed that, indeed, these tracers 
suggest that the density distribution behaves as a power 
law of slope ~ —1, as hypothesized, over at least an order of 
magnitude in density, lending support to our picture of the 
SNRs being largely in their Sedov phases. 

If indeed the end of the visibility of an SNR in the MCs 
is determined solely by its transition to the radiative phase, 
we can, in principle, determine at every point in the Clouds, 
the "visibility time" of SNRs, i.e., the time during which a 
SNR would be visible, if it were there. This variable, (often 
called the "control time" in SN surveys) is an essential input 
to any SN rate calculation. As in Eq. [5] for the transition 
radius, we can derive the dependence of the transition time, 
tmax, on the explosion energy and the ambient density, 

f nu pt 2 - 2 '!)^ 11 -^) -(7-2,0/(11-60 , . 

Again, for a cooling- function power-law dependence on tem- 
perature with index of e = —1/2 to —3/2, the dependence of 
the visibility time on the density is in a limited range, from 
p" 4 / 7 to p- 1 ' 2 . 

To calculate an accurate SN rate, one needs not only 
the dependence of the visibility time on the variables, but 
also the constant of proportionality. For example, a factor 2 
increase in visibility time translates directly to a decrease by 
the same factor of the SN rate. Unfortunately, it is impossi- 
ble to know, except crudely, what is the absolute transition 
time at some point in the MC. The identification of the 
transition to the non-adiabatic phase with the point in time 
when the age equals the cooling time is merely an order-of- 
magnitude device. The proportionality constant in Eq.[7]will 
depend, among other things, on the approximate power-law 
slope e chosen to replace the true cooling function, the nor- 
malization of the cooling function, which is a function of the 
metallicity, the Mach number of the shock, and its adiabatic 
index. Reasonable variations of those parameters alone could 
already change the proportionality constant by an order of 
magnitude. A fiducial value of 10 51 erg is often assumed for 
Eq, the initial kinetic energy of a SN explosion. However, 
there are few SNe or SNRs that have sufficient data and 
are at accurately known distances such that Eo can be reli- 
ably determined. The distribution P(Eo) is also uncertain, 
although energy budget arguments can be made to argue 
that extreme deviations from the fiducial value should be 
rare in both CC and Type la events. Because the depen- 
dence of t max on Eo is mild, roughly -Eg 7 ' 4 or weaker, P{Eo) 
should not be a major source of variations in the mean visi- 
bility time, but the precise magnitude of these variations is 
difficult to quantify. We conclude that, based on the above 
analytic treatment, the mean transition time to the radiative 
stage and hence the mean visibility time of SNRs is known 
only to order-of-magnitude accuracy, (see, however, £16.81 'l 

It would be hoped that the relation in Eq. [7] could be 
calibrated empirically, e.g., if we had an independent age es- 
timate for a SNR of measured size, known to be in its Sedov 
phase. However, there is currently no way to obtain reliable 
age estimates for the relatively large and old SNRs consid- 
ered here (although some constraints can be derived from 



SNR-pulsar associations, see iGaensler fc Johnston! 1 19951 . 
and references therein - we revisit this issue in jj6]) . The dis- 
covery of light echoes associated with SNRs J0509. 0—6844 
(N103B), J0509.5-6731 ( B0509-67.5) and J0519.6-6902 
(B0519-69.0) in the LMC (|Rest et al.l(2005l ) has opened in- 
teresting possibilities for a few objects, but they are all too 
small and young to serve as calibrators. To proceed with our 
calculation of the SN rate in the MCs, we will therefore sim- 
ply leave f ma x,o, the transition time out of the Sedov phase 
for a SNR at a location having the mean ambient gas den- 
sity of the Clouds, as an adjustable parameter. At different 
locations in the MCs, we will scale the visibility time accord- 
ing to the density, as indicated by various tracers at those 
locations, according to 

f — + „ n- A / 7 to - 1 / 2 

^max — fTnax,0 p ■ V J 

This will give us the relative visibility time at each loca- 
tion, which will be useful in i|6] below, when comparing SNR 
numbers and star fomation histories in indvidual MC cells. 
The proportionality constant, i ma x,o will be obtained by re- 
quiring that the CC-SN yield (i.e., the number of CC-SNe 
per unit stellar mass formed) corresponds to the theoreti- 
cal expectation that most massive stars undergo CC-SN ex- 
plosions. This expectation is sup ported by an indep endent 
estimate of the DTD of CC-SNe (|Maoz et al.ll201(M ). 



4 SNR VISIBILITY TIMES FROM THREE 
ESTIMATES OF THE GAS DENSITY 
DISTRIBUTION IN THE MAGELLANIC 
CLOUDS 

In Paper I, we studied three different tracers of gas den- 
sity in the MCs: the HI column density; the star-formation 
rate (SFR) based on resolved stellar populations; and the 
Ha emission-line luminosity. We used these tracers to test 
successfully the hypothesis that the flat distribution of SNR 
sizes results from a transition out of the Sedov phase, which 
is, in turn, determined by the ambient density, combined 
with a density distribution that is approximately a power 
law of index —1. We now use the same three tracers of the 
gas density to scale the SNR visibility time. 



4.1 21 cm emission-line- based HI column density 

The surface brigh tness of HI 2 1 cm line e mission in the MCs 
from t he maps of iKim et al.l (|2003l ) and IStanimirovic et al.l 
|l999f ). analysed in Paper I, is optically thin, and hence 
directly proportional to the HI column density. Since 
the LMC posesses a fairly face-on (inclination i ~ 35°, 
Ivan der Marel fc Cionil l200ll ). well-ordered HI disk, the col- 
umn density should, in turn, be roughly proportional to the 
volume density p. Assuming that the mean volume density p 
in a cell is proportional to its mean column density Nh, we 
can scale the visibility time for each cell containing a SNR 
according to Eq. [8] with t ma x,o assigned to cells having the 
mean Nh of the clouds, 1.5 x 10 21 cm -2 per cell. 
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4.2 Schmidt Law plus star formation rates from 
resolved stellar populations 

As a second way to estimate the local gas densities, 
in Paper I we used the recen t SFRs in each c ell o f 
the lHarris fc Zaritskvl (|2004 ) and lHarris fc Zaritskvl (|2009l ') 
maps. Here, we use this tracer of gas density also to ob- 
tain the scaling of the visibility times of SNRs in the MCs. 
We translate th e 35-Myr-averag ed S FR to a mass colu mn, 
using an inverse ISchmidtl (| 1959T ) law. IKennicuttl (1 19981 ) up- 
dated the Schmidt law, relating star formation rate surface 
density, Esfr, to gas mass column Eg as , as 

/ y, \ 1.4±0.15 

Ssfr = (2.5±0.7) x 1(T 4 - gas M Q yr^kpc" 2 . 

\M e pc~ 2 J 

(9) 

Recent measurements and discussion s of the Schmidt - 
Kennicutt relation can be found, e.g., in iBigiel et al] (|200Sf ). 
and references therein. As emphasized by these studies, the 
Schmidt law has a threshold at some mass c olumn, be- 
low w hich the star-formation rate falls steeply. IKennicuttl 
1 19891 ). for example, finds a threshold at a mass column 
of E gas = 3 — 4 Mqpc -2 , but threshold values several 
times higher have also been reported. From Eq. [5] with 
Egas = 3.5 Mqpc" 2 , the threshold mass column corresponds 
to a SFR of 0.0014 M yr^kpc -2 , or 4.3 x 10~ 5 M yr" 1 
per 12' x 12' cell in the lHarris fc Zaritskvl < |2009l > maps of 
the LMC. For the 24' x 24' cells in the LMC the threshold 
is of course 4 times higher, and for the 12' x 12' cells in the 
SMC, which is 20% more distant, it is 1.44 times higher. We 
therefore use the inverse Schmidt law with an exponent of 
1/1.4 = 0.714 to obtain the density in each cell, down to 
this level of SFR. Below this SFR level, we assign to the cell 
a constant density corresponding to the threshold level. The 
density, in turn, is translated to a visibility time according 
to Eq. [5] with t max ,o assigned to cells having the mean SFR 
of the clouds, 3.3 x 10~ 4 M© yr _1 per cell. 

4.3 Schmidt Law plus star formation rates from 
Ha emission 

Rather than measuring SFR via the resolved stellar popu- 
lations of the MCs, we can trace it by means of Ha emis- 
sion, which then gives us a third tracer of gas density. Ha is 
principally powered by photoionisation from O-type stars, 
whose numbers are proportional to the SFR. The Schmidt 
law then, again, connects SFR and gas mass column. In Pa- 
per I, we analyzed the continuum-subtracted Ha emission 
maps of the LMC an d SMC from the SHASSA survey of 
iGaustad et all (|200lT ). and showed that both the LMC and 
the SMC have distributions of Ha surface brightness that 
follow power laws of index —1, over 2 orders of magnitude 
in flux for the SMC, and over 3 order of magnitude for the 
LMC. Through the Schmidt law, this implied a power-law 
with index —1 for the gas density distribution, over 2 orders 
of magnitude in density. Here, we use the same Ha maps and 
translate the Ha surface brightness to a visibilit y time via 
the SF R and the Schmidt law. We again use the IKennicuttl 
|l989l ) threshold of 4.3 x 10~ 5 M yr _1 per MC cell. By 
comparing the mean Ha surface brightness and the SFR in 
each cell, we find that this threshold corresponds to ~ lOOdR 
(deciRayleighs) in Ha flux per typical MC cell (this holds 
for all cells, both in the LMC and the SMC, since surface 



brightness is independent of distance for these nearby galax- 
ies). Below this Ha flux, we assign to a cell, in the absence 
of a density indicator, a constant density corresponding to 
the threshold level. 



5 RECONSTRUCTION OF THE SN DELAY 
TIME DISTRIBUTION - METHOD 

With our SNR sample, and our estimate for the visibility 
time of SNRs at each location in the MC in hand, we now 
briefly present a method to combine this information with 
the localized SFHs for each individual cell in the MC, in 
order to recover the SN DTD. A more detailed exposition 
of the method, including demonstrati ons of its performan ce 
on simulated datasets, is presented in lMaoz et al. I l|2010ah . 

The SN rate in a galaxy observed at cosmic time t is 
given by the convolution 

r(t) = f S(t-T)*(r)dT, (10) 
Jo 

where S(t) is the star-formation rate (stellar mass formed 
per unit time), and \&(t) is the DTD (SNe per unit time per 
stellar mass formed). As reviewed is §1, previous attempts 
to recover the DTD have used rates r(t) measured in surveys 
of galaxies at different redshifts (i.e., different cosmic times), 
compared to cosmic star-formation histories, whether in field 
surveys or galaxy cluster surveys. An alternative approach 
has been to look at the SN rates per unit stellar mass in 
galaxies of particular types (star-forming, quiescent, etc.), 
and to attempt to assign to each type a "formation age" or 
some generic, simple, star-formation history. A shortcoming 
of all these approaches is that they involve averaging over 
the galaxy population (i.e., all the SNe are assumed to come 
from the entire host population considered), or over time 
(i.e., the detailed history of all galaxies of a certain type is 
represented by a single "age" or simplified history). As a 
result, all of these approaches involve loss of information, 
and potential systematic errors (e.g., due to unrepresenta- 
tive simplified histories). 

An alternative way to recover the DTD is by inverting a 
linear, discretized version of Eq. 1101 where the detailed his- 
tory of every individual galaxy or galaxy subunit is taken 
into account. Suppose the star-formation histories of the 
i — 1,2, N galaxies or galaxy subunits monitored as part 
of a SN survey are known (e.g., based on reconstruction of 
their stellar populations), with a temporal resolution that 
permits binning the stellar mass formed in each galaxy into 
j — 1,2,.., K discrete time bins, where increasing j corre- 
sponds to increasing lookback time. The time bins need not 
necessarily be equal, and generally will not be, since the tem- 
poral resolution of the star-formation-history recontruction 
degrades with increasing lookback time. For the ith galaxy 
in the survey, the stellar mass formed in the jth time bin is 
mij . The mean of the DTD over the jth bin (corresponding 
to a delay range equal to the lookback-time range of the jth 
bin in the star- formation history) is Then the integral 
in Eq. [TU]can be approximated as 

K 

-<Y,"' Al - ( n ) 

where n, the SN rate in a given galaxy, is measured at a 
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particular cosmic time (e.g., corresponding to the redshift 
of the particular SN survey). Given a survey of N galaxies, 
each with an observed SN rate n and a known binned star- 
formation history rriij, one could, in principle, invert this 
set of linear equations and recover the best-fit parameters 
describing the binned DTD: * = * 2 , 

In practice, on human timescales SNe are rare events in 
agiven galaxy, i.e., n <C 1 yr _1 . Supernova surveys therefore 
monitor many galaxies, and record the number of SNe dis- 
covered in every galaxy. For a given model DTD, Vfri, .., *&k, 
the ith galaxy will have an expected number of SNe 

Aj = nu, (12) 

where ti is the effective visibility time during which a SN 
would have been visible (given the actual on-target moni- 
toring time, the distance to the galaxy, the flux limits of the 
survey, and the detection efficiency). Since A; <C 1, the num- 
ber of SNe observed in the ith galaxy, m, obeys a Poisson 
probability distribution with expectation value Xi, 

P(n i \\ i )=e- Xi X? i /r H l, (13) 

where m is zero for most of the galaxies, 1 for some of the 
galaxies, and more than 1 for very few galaxies. 

Considering a set of model DTDs, the likelihood of a 
particular DTD, given the set of measurements ni,...,njv, 
is 

L = l[P(ni\Xi). (14) 

i=l 

More conveniently, the log of the likelihood is 

JV JV JV 

lnL = ^lnP(ni|Ai) = - ^ Xi + ^ m(A™V^!), (15) 

i— 1 i— 1 i — 1 

where obviously only galaxies hosting SNe contribute to the 
second term. The best-fitting model can be found by scan- 
ning the parameter space of the vector for the value that 
maximizes the log-likelihood. This procedure naturally al- 
lows restricting the DTD to have only positive values, as 
physically required (a negative SN rate is meaningless). 

The covariance matrix Cjk of the errors in the best- 
fit parameters can be found by calculating the curvature 
matrix, 

_1 G> 2 lnL A 9[lnP(n i |A I )] d[ln P(n, | A,)] 
J i=l J 

and inverting it, 

[C] = [a]' 1 . (17) 

Because the values of the DTD are constrained to be posi- 
tive, if the maximum- likelihood of a DTD component 9j is 
close to zero, the square-root of its variance, yj Cjj , will not 
represent well its la uncertainty range. An alternative, more 
reliable, procedure is to perform a Monte-Carlo simulation 
in which many mock surveys are produced, each having the 
same galaxies, star-formation histories, and visibility times 
as the real survey, and having expectation values Xi based 
on the best-fit DTD, but with the number of SNe in ev- 
ery galaxy, m, drawn from a Poisson distribution according 
to Xi. The maximum-likelihood DTD, \E>, is found for ev- 
ery realization. From the distribution of the values of every 



component, tyj, over all the realizations, one can estimate 
the range encompassing, say, 134% of the cases. Another 
advantage of the Monte-Carlo simulations is that they allow 
gauging the effects of additional sources of error, such as 
uncertainties in the SFHs. 

The above approach for recovering the DTD has sev- 
eral advantages over previous methods. First, all the known 
information in the survey is included in the analysis in a sta- 
tistically rigourous way, including the fact that many (usu- 
ally most) of the galaxies, or galaxy subunits, did not host 
any SNe. Furthermore, the calculation is easily generalized 
to cases where the galaxies are not all at the same distances 
(e.g., combinations of surveys done at different redshifts); 
one simply needs to use the appropriate star-formation his- 
tory bins for every galaxy. In fact, assuming the DTD is 
a universal function (i.e., it is independent of environment, 
metallicity, cosmic time, etc., see §1) it is straightforward to 
combine, in a single analysis, the data from completely dis- 
parate SN surveys, e.g., classical SN surveys, where a large 
volume or a lare sample of galaxies is monitored, with un- 
conventional SN "surveys" such as the one presented here, 
in which the SN rate is measured based on SN remnants in 
small subunits of a single galaxy. We re-emphasize that, in 
the preceding discussion, everything relating to "galaxies" 
is equally applicable to individual subunits of a galaxy, for 
which a SFH is available. 

As described in more detail in iMaoz et all (|2010ah . the 
number and resolution of the time bins used in the analysis 
will naturally depend on the quality of the data; the larger 
the number of observed SNe, -/V to t, the higher the time res- 
olution that can be recovered with reasonable accuracy. For 
a survey with a fixed -/V to t, there will be a tradeoff in the 
analysis between DTD accuracy and resolution. 

6 THE DTD IN THE MAGELLANIC CLOUDS 
6.1 General considerations 

We now apply the IMaoz et all (|2010al l DTD recovery 
method, described in the previous section, to the 77 SNRs 
in the Clouds, treating the sample as a SN survey with visi- 
bility times given by Eq.[8l co nducted over the 1836 gal actic 
subunits or "cells" define d bv lHarris fc Zaritskvl (|2004l ) and 
lHarris fc Zaritskvl (|2009l 'l. each with its reconstructed SFH, 
based on the resolved stellar populations. We recall that 
changing the mean visibility time, t m ax,o, which still needs 
to be adjusted, causes the amplitude of the DTD to scale 
identically across all time bins, as i^ax o- 

One point of concern in the derivation of the DTD could 
be that, because of the random stellar velocities in each 
galaxy, the stars currently in a given cell which contributed 
to the SFH are not the same stars that produced the SNe 
observed to have exploded in that cell, many years after the 
formation of their progenitors. While this is true, it is in- 
consequential to a correct determination of the global DTD, 
because the same spatial diffusion affects both the progen- 
itor population and those stars among it that eventually 
explode. To see this, consider, as an example, a grid of 3 x 3 
adjacent cells in the LMC. Suppose, as a toy example, that 
500 Myr ago there was a short burst of star formation in the 
central cell, forming a stellar mass M, and no activity in the 
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other cells. Suppose, further, that the SN DTD is such that 
the stellar population formed in the burst leads, 500 Myr 
later, to nine SNe over the past 20 kyr, which are therefore 
detected as SNRs in the galaxy, i.e., a ratio of 9/M SNe per 
unit stellar mass formed 500 Myr ago. Finally, suppose that 
the stellar diffusion timescale in the galaxy is such that, over 
the 500 Myr, the progenitors of the 9 SNe, before explod- 
ing, have drifted out of the central cell in which they were 
formed, and there is now, on average, one SN in each cell. 
However, the entire stellar population of the burst will have 
diffused in the same way, and therefore each cell will have 
1 /9 of the 500 Myr old population that was originally in the 
central cell. When we compare SN numbers to the 500 Myr- 
old stellar mass present today in each cell, we will see 1 SN 
per M/9 of stellar mass formed, and therefore we will still 
deduce the correct ratio of 9/M SNe per unit stellar mass 
formed. This argument holds no matter what are the look- 
back times, the diffusion timescales, or the complex SFH in 
each cell. 

A peculiarity of our SN survey is that, for most of the 
remnants, we have little or no information on the type of the 
SN that exploded. Some young SNRs can be classified us- 
ing different methods: SNR 1987A was obviously a CC-SN, 
SNR J0509. 5-6731 (B0509-6 7 5) was a SN la b oth from the 
spectroscopy of its light echo dRest et al. 20081) and t he X- 
ray emission from the SN ejecta ( Badenes et al.l [20081 ). and 
several other objects have classifications that range from 
the very secure (those harbouring compact objects are al- 
most certainly CC, see below) to the reasonable (young ob- 
jects with ejecta-dominated X-ray emission). However, these 
methods usually cannot be applied to the old SNRs in the 
Sedov phase that form the bulk of our obje cts (for a more 
detaile d discussion on SNR typing, see § 3 in lBadenes et al] 
l2009al ). The SNR sample therefore constitutes a mix of dif- 
ferent SN types. However, this is easily addressed within our 
formalism. As already noted above (see[T] CC SNe will gen- 
erally explode within tec ~ 30 — 40 Myr of star formation, 
and SNe la will explode after this time. We will therefore 
choose our shortest time bin as < r < 35 Myr, and inter- 
pret the DTD amplitude in this bin as the CC SN signal, 
while the signal in later bins is due to the SNe la. 

A systematic error arises through the uncertainty in the 
visibility time, due to the the cooling function dependence 
on temperature (Eqns. I7I8[) . and due to the three different 
tracers of density that we consider - HI column density, re- 
solved SFR, and Ha luminosity (3J). We probe the effect of 
this systematic uncertainty, by recovering the best-fit DTD 
for all combinations of the density tracers with the two ex- 
treme density exponents of Eq. [8] —4/7 and —1/2. We find 
that the choice of exponent affects the derived DTD ampli- 
tudes only slightly, by < 3% for the HI and Ha tracers, and 
by 10% for the resolved star formation rate tracer. We will 
henceforth cite results only for an exponent of —1/2. 

Figure [T]and Table 2 show the Magellanic Cloud DTD 
we recover from the SNR sample, binning the DTD into the 
following three intervals: r < 35 Myr (CC SNe), 35 Myr < 
r < 330 Myr ("prompt" SNe la) and 330 Myr < r < 14 Gyr 
("delayed" SNe la). The small number of SNe in our "sur- 
vey" precludes the possibility of finer binning and a detailed 
recovery of the DTD shape. Nevertheless, we will see that 
some interesting results emerge even with this coarse bin- 
ning scheme. 



6.2 Core-collapse SNe 

For each of the density tracers, we renormalise the vis- 
ibility time at the mean density, i ma x,o, such that the 
rate in the first DTD bin, which traces CC-SNe, is = 
2.86 SNe yr -1 (10 10 M a ) -1 . Multiplied by the width of the 
bin, 35 Myr, this then gives a time-integrated CC-SN yield 
(i.e., the number of CC-SNe per unit stellar mass formed) 
of Ncc/M = 0.01 M _1 . This is the value expected if all 
stars above 8Mq explode as CC SNe, 



J™°{dN/dm)dm 



Ncc . 

M 0.7 f™°m(dN/dm)dm 



= 0.01 M 



© 



(18) 



for a "diet Salpeter" IMF (|Bell et al.ll2003l ), where dN/dm oc 
m~ 2,35 , m is stellar mass in units of Mq, and the factor 
0.7 in the denominator accounts for the reduced number of 
low-mass stars in a realistic IMF, compared to the original 
ISalpeteij (| 19551 ) IMF (see 

As listed in Table 2 for the three density tracers, the 
renormalized visibility time at a location in the MCs hav- 
ing the mean density is: 22.5 kyr (resolved SFR); 13.9 kyr 
(Ha); and 13.3 kyr (HI). Depending on the local value of 
the density, individual SNRs in the MCs could obviously be 
much older - at a density 10 times lower than the mean, for 
instance, the SNR lifetime would be a factor 3 to 4 longer, 
depending on the value of the exponent in Eq. [S] Our SNR 
lifetimes are therefore in agreement with the maximum ages 
of S NR shells constra ined by SNR-pulsar associations (~ 60 
kvr. lFrail et al.lll994h . 

Our resulting DTD will be incorrectly normalized if, 
in reality, the mass border between CC-SNe and SNe la 
is not at 8A/q, or if a large fraction of massive stars end 
their lives without a SN explosion that l eaves a remnant 
(e.g.. through direct collapse to a black hole; Kochanek et al.l 
[20081) . However, a direct derivation bv lMaoz et al.l ( 2010a ) 
of the CC-SN DTD from the LOSS-SDSS sample supports 
the conclusion that the CC-SN yield is at the level assumed 
here. 

The errors we cite in Fig. \T\ and Table 2 give the 68% 
probability range of each component, based on Monte-Carlo 
simulations, as described above. In these simulations^ to ao 



count also for the uncertainties giv en by lHarris fc Zaritskvl 
|2009l ) and lHarris fc Zaritskvl (|2004l ) for their SFHs, we have 
done the following. For each MC cell, we add in quadrature 
the la errors in SFH among the fine time bins that consti- 
tute a single co-added time bin in our analysis. Since there is 
covariance among the SFH errors in adjoining bins, this ad- 
dition constitutes a conservative overestimate of the true er- 
rors. Then, in every mock survey in the simulation, we draw, 
for each cell in every time bin, a SFR from an asymmet- 
ric Ga ussi an distribution, cen t ered o n the lHarris fc Zaritskvl 
l|2009h and lHarris fc Zaritskvl (|2004l ) best-fit value, with pos- 
itive and negative standard deviations according to the bin- 
added negative and positive errors above. After calculation 
of the expectation value of the rate for the cell, according to 
Eq. 1121 the "observed" number of SNRs in the cell is drawn 
from the appropriate Poisson distribution. An exception to 
this procedure is in the SM C, where we have avoided us- 
ing the la positive errors of lHarris fc Zaritskvl (|2004l ). as 
they are two orders of m agnitude larger than those given by 
lHarris fc Zaritskvl (|2009j ) for the LMC, despite the similarity 
between the two galaxies, their data, and their analysis. For 
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Figure 1. Best-fit delay-time distributions for the Magellanic Cloud SNR sample. Each panel in the top row displays the results obtained 
with a different density tracer: HI (left), Ho (centre), and SFR (right). The panels in the bottom row show the DTDs obtained with 
"hybrid" density tracers, using HI at low densities and either Ho (left) or SFR (right) at high densities. Error bars show 68% confidence 
intervals, from Monte-Carlo simulations that account both for Poisson statistics and for the uncertainties in the SFR of each cell. The 
best fits in the most-delayed bins are always zero, and the horizontal line gives the 95% confidence upper limit on the rate in this time 
bin. To facilitate comparison, the DTD obtained without scaling with density of the visibility time of SNRs is plotted in all panels as 
black lines, with shaded regions for the error bars. All DTDs have been scaled so as to give the same, theoretically expected, CC-SN 
yield in the first bin. 



the SMC, we therefore use the negative errors in the SFH 
to represent the positive errors as well. Finally, in the DTD 
reconstruction of each mock survey, the Harris & Zaritsky 
(2004, 2009) best-fit value is assumed for the SFH (i.e., a 
value generally different than what was used to generate the 
mock survey). 

We note that, although we have renormalized t m ax,o, 
so as to always obtain *i = 2.86 SNe yr -1 (10 10 A/ Q ) _1 , *i 
is still a parameter in the DTD reconstruction, and its de- 
rived error indicates the confidence level at which a CC-SN 
component is detected in the DTD. For the different density 
tracers, the ^\ component of the DTD is detected in our re- 
construction at the 5 — 7a level. In other words, there is no 
doubt that some of the SNRs in the MCs are significantly 
associated with stellar populations that are young enough 
to produce CC-SNe. 



6.3 Prompt SNe la 

The strong DTD signal recovered in the CC bin effectively 
"uses up" much of the SN data in the survey, leaving a 
weaker signal in the SN la bins. Nevertheless, there is a sig- 
nal in the "prompt" SN la bin, * 2 , at 35 Myr< r < 330 Myr, 
at the 1.7a level for the resolved star formation tracer, 2.3er 



for the Ha tracer, and 3a when using HI column as the den- 
sity tracer. The rate ^2 itself is also correspondingly higher 
or lower when using each of these tracers, as shown in Ta- 
ble 2. Here, by u Xa" we mean the number of —34% uncer- 
tainty intervals above zero. Because of the non-Gaussian na- 
ture of the errors, the significance of the detection is higher 
that it would be for the same number of standard devia- 
tions in a Gaussian. Specifically, we find from Monte-Carlo 
simulations, in which we input a value of $2 = 0, that re- 
covered values as high as the best-fit values are obtained in 
fewer that 1% of mock surveys, using any of the tracers. In 
other words, a nonzero prompt SN la component is detected 
at > 99% confidence. Using additional simulations, we have 
determined, for each tracer, the input ^2 rate that gives re- 
covered values greater of equal to the best-fit values in 
5% of the realizations. This input ^2 rate is the 95% con- 
fidence lower limit on $2, which is also listed in Table 2. 
Despite the clear detection of a prompt SN la component, 
there is factor 4 systematic uncertainty in the best-fit am- 
plitude of this component. This stems from the systematic 
uncertainty in estimating the local density, which in turn 
leads to an uncertainty in the scaling of the visibility time. 

The $2 rate, multiplied by the length of the second 
time bin (295 Myr), indicates a time-integrated produc- 
tion of prompt SNe la of: Ni & /M = 0.0027 ± 0.0014 M _1 
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Table 2. Mean SNR Visibility times, SN rates, and reconstructed DTDs from the MC SNR sample. 



Density 
Tracer 


Mean SNR 
Visibility 
[kyr] 


SN Rate 
[1(T 3 SN yr^ 1 ] 


[SNuM] 


0-35 Myr 


*i [SNe yr" 
35-330 Myr 


- 1 (io 10 Mg)- 1 ] 

*2 

(95% low lim.) 


*3 

0.33-14 Gyr 


HI 

SFR (Schmidt law) 
Ha (Schmidt law) 
No Scaling 
HI+Ha 
HI+SFR 


13.3 
22.5 
13.9 
13.4 
12.6 
15.9 


4.1 ±0.9 
2.4 ±0.4 
3.3 ±0.6 
5.8 ±0.1 
3.8 ±0.8 
3.1 ±0.6 


3.0 ±0.7 
1.7 ±0.3 
2.3 ±0.4 

3.5 ±0.7 

2.6 ±0.5 

2.1 ±0.4 


2 -° b ~0.50 

2.86ii» 
2.861°;" 


0.389 ±0.131 
0.091 ±0.057 
0.194 ±0.086 
0.316 ±0.125 
0.260 ±0.100 
0.165 ±0.080 


(> 0.171) 
(> 0.016) 
(> 0.060) 
(> 0.143) 
(> 0.104) 
(> 0.058) 


< 0.0014 

< 0.0014 

< 0.0016 

< 0.0024 

< 0.0014 

< 0.0013 


SN rates in units 


of SN yr" 1 


are for the LMC 


and SMC 


together. 


SN rates per unit stellar mass 


are in SNuM 



[SNe (100 yr) - (10 Mq)~ ], after converting formed mass to present mass, by assuming 30% mass loss. Quoted num- 
bers are: best-fit values and 68% confidence intervals for the rates, \E' i and >I/2; 95% confidence lower limits for >l/2; and 95% 
confidence upper limits for ^3. 



a The rate in this bin was renormalized, but the errors reflect its significance, see text for details. 



(resolved SFR); Nu/M = 0.0057 ± 0.0026 Mo" 1 (Ha); 
or JVia/M = 0.0114 ± 0.0039 M Q _1 (HI). We note that 
Ni a/M, above, is essentially the " B parameter" di s cusse d 
by IScannapieco fc Bildstenl |2005l ). iMannucci et alj (|200rJ ). 
and ISullivan et al.l | 2006T k i.e.. the ratio of SN la rate to 
SFR i n stro ngly star-forming galaxies. As summarised in 
iMaozl |2008l ) , the values of B generally found by SN surveys 
are in the range B = 0.0 01 to 0.003 Mq -1 , an exception be- 
ing [Sulliva^et^al] (|2006l ). who obtain a value several times 
lower, (3.9±0.7) x IO^Mq" 1 . The values we find for Nu,/M 
are broadly consistent with those estimates, but tend to be 
on the high side, particularly when we use HI column as a 
density tracer in the MCs. 

The time-integrated ratio of CC-SNe to SNe la (or 
equivalently, the ratio of rates of CC and type-la SNe in 
a population that has had a constant star-formation rate 
for a long time, such that both rates have reached a steady 
state) is in the range of 9 : 1 to 1 : 1, when we consider both 
the statistical and the systematic uncertainties. 



6.4 Delayed SNe la 

The DTD amplitude in the third bin has a best-fit value of 
zero in all of our reconstructions, but an uncertainty range 
that includes the typical SN la rates that have been mea- 
sured in old, quiescent populations. The uncertainty is large 
because few SNe are expected in this bin in our small sam- 
ple, given the dominance of young stellar populations in the 
MCs. Specifically, from our Monte-Carlo simulations, the 
95% confidence upper limit on the delayed SN la rate in 
the DTD is * 3 < 0.0013 to 0.0016 SNe yr- 1 (10 10 M Q )- 1 
(depending on the density tracer). In other words, such an 
input value in the simulations for this component of the 
DTD results in a best-fit reconstructed value of zero in 5% 
of cases. For comparison, the typical values found for the 
11 A parameter" , the SN la rate per unit mass in an old pop- 
ulation that has no ongoing star formation, have been in 
the range A ~ (2 - 10) x 10~ 4 SNe yr _1 (10 10 M s ) _1 (see 
compilation in lMaozll200g| ). However, ^3 is the rate per unit 
mass formed. The published rates in old populations are 
a per unit of existing stellar mass, i.e., in stars and stel- 
lar remnants. In the course of ~ 10 Gyr, a stellar popu- 



lation (with a "diet Salpeter" IMF that we are assuming 
throughout) will return w 50% of its mass to the interstellar 
mediu m via SN explosions and m ass loss during stellar evo- 
lution (|Bruzual fc Charlotl Il993l ) . Thus, for comparison to 
^3, A needs to be multiplied by 0.5, giving ~ (1 — 5) x 10 -4 
SNe yr -1 per (1O 1O M ) -1 formed. This range fits comfort- 
ably, by at least a factor of 2, below the upper limit we 
have derived here based on the SNRs in the MCs. An ap- 
plication of our DTD recovery method to the LOSS SN sur- 
vey, which does hav e a significant old galaxy population, by 
iMaoz et al.l (|2010al ). gives a highly significant measurement 
of a delayed SN la component, * 3 = 2.6±g;| x 10~ 4 SNe 
yr _1 (10 M©) . This value is consistent both with previ- 
ous determinations of the A parameter and with the upper 
limit on ^3 found here. 



6.5 Relative fraction of prompt and delayed SNe 
la 

Bins 2 and 3 of the DTD thus, respectively, require the exis- 
tence of prompt SNe la, and allow (but do not require) the 
existence also of a delayed component. We note that, if we 
choose to ignore the density distribution, and assume a con- 
stant visibility time throughout the MCs, we obtain results 
that are intermediate to those using Ha and HI, in terms 
of visibility time and prompt-SN-Ia significance and ampli- 
tude, and hence also in terms of SN la yield and CC-SN to 
SN la time-integrated ratio. 

We can also use the ^2 measurements and the upper 
limits on *l/ 3 to gauge the relative fraction of the prompt 
and delayed SNe among the SN la population. Using the 
resolved-SFR-based values of ^2 and ^3 in Table 2, which 
have the lowest ratio, we can obtain a lower limit on the 
time-integrated SN la yields in the two components, 



*2At2 

^3 At. 



3 



>1, 



(19) 



wher At2 and At 3 are the sizes of time-bins 2 and 3. Using 
the other tracers gives lower limits on this ratio that are 
larger. In any case, a significant fraction of SNe la, possibly 
a majority, are prompt. 
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6.6 Hybrid density tracers 

Among the three density tracers that we use, the HI column 
density is physically closest to the actual volume density 
we are interested in. Furthermore, it can probe and repre- 
sent the lower-density regions of the MCs, where the two 
other tracers, which rely on an inverse Schmidt law, break 
down because of the mass column threshold in the SFR. In 
those low-density regions, we have knowingly assigned a con- 
stant, artificially high, density which systematically short- 
ens the visibility time. Finally, the Schmidt law has a large 
intrinsic scatter, inducing noise in the assignment of visi- 
bility times. If we focus on the Hi-tracer results, we find a 
highly significant prompt- SN-Ia component, with an ampli- 
tude somewhat higher than found by field SN surveys, and a 
time-integrated ratio of CC-SNe to SNe la of approximately 
1:1. Interestingly, by looking at this ratio for the youngest 
remnants in the LMC, which are the ones that can be most 
securely classified, one would have guessed as much - four 
of the eight sm allest SNRs in the L MC were definitely or 
likely type la's (|Badenes et alj|2009al v 

On the other hand, HI as a density tracer also has its 
shortcomings. For example, the cell in the centre of the 30 
Doradus region of the LMC has a high SFR, as evidenced by 
its resolved mean SFR over the last 35 Myr, which is 5 times 
higher than the mean of the MCs, and by its Ha luminos- 
ity, which is 30 times higher than the mean. The HI column 
densi ty, however, is only twice as high as the mean for the 
MCs. lHarris fc Zaritskvl (|2009T ) have already noted the far- 
from-perfect spatial correlation amongst different tracers of 
star formation. In the case of 30 Doradus, the HI column is 
likely underestimating the true total gas density, because a 
fraction of the neutral hydrogen has been photoionised by 
the massive stars in this region. An underestimate of gas 
density in such regions will lead to an overestimate of the 
visibility time. This, in turn, will lower the DTD rate 
of CC-SNe that are associated with such regions. Renor- 
malising the entire DTD such that = 2.86 will then lift 
also the prompt-SN-Ia component, ^2, causing its rate to be 
overestimated. At some level, this is likely occurring when 
we use HI column as a density tracer. Furthermore, we have 
discussed in Paper I the evidence that HI may be misrepre- 
senting the density in the SMC because of projection effects. 
Thus, this tracer is also not foolproof. 

To deal with these problems, we have combined two 
different tracers, using each in the regime where it is likely 
to be most reliable. We have rederived the visibility time 
in each cell, using HI column as the density tracer in cells 
where the column density is lower than the mean (about 
2/3 of the cells), and using SFR and the inverse Schmidt 
law, as before, in the rest of the cells, where HI is above 
the mean. The SFR is traced either with the resolved stellar 
populations, or using Ha. 

The best-fit DTD recovered with these visibility times 
(see Table 2), not surprisingly, is intermediate to the ones 
found with HI and Ha or the resolved SFR when used alone 
as tracers. The prompt SN la component is detected at a 
level of 2.6a above zero (using HI + Ha) or 2.1a above zero 
(using HI + resolved SFR). As before, from the Monte Carlo 
simulations, the significance of these detections is > 99%. 
The results do not depend strongly on the adopted border 
for the use of the two tracers. For example, setting the border 



between use of HI and resolved SFR at 1/2 of the mean HI 
column or at twice the mean HI column (instead of at the 
mean), changes the best-fit value of ^2 by —27% and +20%, 
respectively. We consider these hybrid-tracer-based DTDs to 
be our most reliable results. 

6.7 Sensitivity to temporal binning and leaks 

A concern with all of our results is that they may be sensitive 
to the temporal border between CC-SNe and SNe la, which 
we have set at 35 Myr. For example, if some CC-SNe ex- 
plode, say, between 35 Myr and 45 Myr, some of the CC-SN 
signal will be attributed incorrectly to SNe la. As a result, we 
will overestimate the number of prompt SNe la. To test this 
possibility, we have performed the following experiment. We 
have selected 11 LMC SNRs from Paper I that are almost 
certainly remnants of CC explosions, based on one or more 
of the following criteria: detection of pulsars within them, 
X-ray analysis of the ejecta showing clear CC-SN products, 
or (in the case of SN1987A) direct historical evidence. With 
this small sample of known CC-SNe, which we list in Table 3, 
we have re-derived the DTD. Naturally, the amplitude of the 
shortest, CC-SN, bin of the DTD will be understimated, be- 
cause we have selected only a fraction of the CC remnants in 
the LMC. Nevertheless, we can test if any signal "leaks" into 
the later bins, which we have been associating with SNe la, 
either because of the time border issue discussed above, or 
becau se of some other systematics of the lHarris fc Zaritskvl 
|2009l ) SFH reconstruction. We find no evidence of "leakage" 
outside the < 35 Myr bin. For example, using the HI+Ha 
hybrid density tracer, the best-fit results are, in the first 
(earliest) time bin of the DTD (with no renormalisation of 
this rate), * t = 0.61 ± 0.26 SNe yr -1 (10 10 M Q ) _:l , and zero 
in the second and third bins. This is reassuring, but because 
of the small size of this subsample, the 95% confidence upper 
limit on the rates in the second and third bins are not low 
enough to rule out, using this test, a leak of the CC signal, 
having the amplitude of the second and third bins in the full 
sample. 

As another test, we have rederived the full-sample 
DTDs with the various tracers, but enlarging the first time 
bin from 0-35 Myr to 0-80 Myr (this step size is dictated 
by the time bins of the Harris and Zaritsky SFH recon- 
structions). Not surprisingly, the best-fit amplitude of \&2 
decreases by factors of 2-3, comparable to the increase in 
the length of the first time bin. Unfortunately, the relative 
uncertainty in ^2 also becomes very large, precluding the 
possibility of saying with any confidence whether or not a 
SN la component at delays of 80-330 Myr exists. However, 
it is highly unlikely that, at 80 Myr delays, there are still 
contributions to the DTD from CC-SNe, as this delay cor- 
responds to the lifetime of stars with initial masses below 
6Mq. Thus, this exercise suggests that much of the signal 
in the 35-330 Myr bin of the DTD could be due to SNe la 
that explode within < 80 Myr. This result joins other ev- 
idence for prompt SNe la with e xtremely short delays, of 
< 100 Myr (jMannucci et al]|200d . ; Delia Valle et al. 2005). 
The obstacle to obtaining a clearer answer, however, is the 
small number of SNRs. One needs to recall that, among the 
77 SNRs in our sample, at most half and possibly only about 
10% are "driving" the SN la signals in the second and third 
DTD bins, and hence the large uncertainty. A larger sample 
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Table 3. Bona fide CC-SNRs in the LMC 



SNR 




Alternate 


CC Classification 






Name 


Criteria and Reference 0- 


J0453.6- 


-6829 


B0453-685 


NS; Gaensler et al. (20031 


J0505.9- 


■6802 


N23 


NS; Hughes et al. ("20061 


J0525.1- 


■6938 


N132D 


EJ; Borkowski et al. (20071 


J0525.4- 


■6559 


N49B 


EJ; Park et al. (2003b1 


J0526.0- 


■6604 


N49 


EJ b ; Park et al. (2003a1 


J0531.9- 


■7100 


N206 


NS; Williams et al. (20051 


J0535.5- 


■6916 


SNR1987A 


historical CC-SN 


J0535.7- 


■6602 


N63A 


EJ; Warren et al. (20031 


J0536.1- 


■6734 


DEM L241 


NS; Bamba et al. (20061 


J0537.8- 


■6910 


N157B 


NS; Chen et al. (20061 


J0540.2- 


-6919 


B0540-693 


NS; Kirshner et al. (19891 



a EJ: Ejecta emission clearly indicative of a CC SN origin (usually, O-rich and Fe-poor); NS: SNR contains either 
a confirmed neutron star or strong evidence for a pulsar wind nebula. 

^ This SNR contains a magnctar, but the NS-SNR connection is disputed (see discussion in lBadenes et al. 2009a). 
Nevertheless, the ejecta composition leaves little doubt about the SN type. 



of SNRs could be obtained by combining SNRs from several 
nearby galaxies, as briefly discussed in §8, below. 



6.8 Comparison of the visibility time to 
theoretical estimates 

An interesting question that we can ask, at this point, is 
how do the observed values of visibility time that we have 
obtained compare to the theoretical expectations for the age 
at the SNR transition to the radiative phase. Blondin et al. 
(1998) found, using ID hydrodynamical simulations and a 
full cooling curve (rather than a power-law approximation) 
that the transition time is well approximated by 1.6 t CO oi for 
a wide range of conditions, where t coo i is their rough estimate 
for the cooling time, as in our Eqns.[2]and[7] assuming a cool- 
ing function represented by an e = —1 power law. Specifi- 
cally, they found that for an ambient hydrogen number den- 
sity of 1 cm" 3 , and for the fiducial explosion energy and 
gas abundances, t coo i m 29 kyr, and hence the theoretical 
transition time for this density is i max ~ 29 x 1.6 = 46 kyr. 

Unfortunately, we cannot directly compare this number 
to the visibility times we have obtained in this work, us- 
ing the various density tracers. Our tracers are all based on 
column densities, rather than densities, under the assump- 
tion that the two are correlated. For example, although we 
know that for the mean HI column density of the clouds, 
1.5 x 10 21 cm -2 , the visibility time is 13.3 kyr (using HI as 
a tracer), we do not know what absolute density this cor- 
responds to. Therefore, neither do we know how to scale 
the visibility time to its value for the fiducial Blondin et al. 
(1998) density of 1 cm" 3 . 

Instead, we will check if the LMC disk thickness needed 
in order to match the numbers is reasonable. Assuming a 
visibility-time dependence on density of i ma x oc p~ 1//2 , our 
observed range of t max = 13.3 to 22.5 kyr implies a mean 
density of n ~ [46/(13.3 to 22. 5)] 2 « 4 to 12 cm" 3 . We note 
that, although Blondin et al. (1998) consider the case of a 
shock expanding into a neutral hydrogen ambient medium, 
the results would change little if the ambient medium con- 
sisted of molecular hydrogen, or a mixture of neutral and 
molecular hydrogen. The 21 cm emission traces only the 



atomic hydrogen, while molecular hydrogen is present in the 
gas disk as well. Norikazu et al. (2001) have used CO as a 
molecular gas tracer to measure a molecular-to-atomic mass 
ratio of 0.2 ; 1 in the low-radial velocity component of the 
gas in the LMC disk, associated with the densest regions. 
However, for low columns or low metallicities, CO may no 
longer be a be a good tracer, because it can be photodisso- 
ciated whereas H2 remains self-shielded (Leroy et al. 2009), 
and the actual molecular ratio is likely higher. Thus, it is 
reasonable to assume that the total gas column density is at 
least « 2 x 10 21 cm~ 2 . Dividing by n then gives a LMC gas 
disk thickness of about 60 to 170 pc or more, smaller than, 
but comparable to the 300 pc-thick Milky Way gas disk. We 
conclude that, while there is too much missing information 
about the LMC geometry (and certainly about that of the 
SMC) for a definitive test of the Blondin et al. (1998) model 
and its application by us to SNR visibility, there appears to 
be rough consistency in the numbers. 



7 THE CURRENT SN RATE 

We have focused, up to this point, on deriving the SN DTD, 
which has permitted us to separate, based on delay times, 
the contributions of CC-SNe, prompt SN la, and delayed 
SNe la. However, the data allow deriving also the more tra- 
ditional current SN rate in the MCs, which is also of inter- 
est. Naturally, the rates we will get will be total rates, for 
all types of SNe, CC and type-la, that produce SNRs. The 
SN rate in the LMC or SMC, or in both galaxies treated as 
a whole, will simply be 



R 



A^SNR-ZVcolls 



(20) 



where Asnr is the total number of remnants, N ce \i s is the 
number of cells, and the sum over the visibility times, ti, is 
over all cells. Alternatively, one can calculate the SN rate 
per unit mass, 



AsNR 
J2i m iti 



(21) 
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where m« is the time-integrated stellar mass formed in each 
cell (i.e., rrii = X^ m ij)- The relative error in these rates 
is given by the sum in quadrature of the relative Poisson 

1/2 

error, N SN ^ , and the relative error of the visibility time. 
Since all the visibility times have been renormalised to give 
*?i = 2.86, the relative visibility-time error is just the rel- 
ative error in vPi. Measurements of Rm in other galaxies 
have generally been normalised not to the formed stellar 
mass, but to the existing stellar mass. For ease of compar- 
ison to those measurements, we will therefore scale down 
all values of rrii by 0.7, which accounts roughly for the 
mass loss of stars during stellar evolution in actively star - 
forming galaxies like the MCs (jBruzual fc Charlotl 1 19931 ). 
(This 30% mass loss is in contrast to the 50% mass loss 
previously considered for old and quiescent stellar popula- 
tions). Rates per unit mass in Table 2 are in units of SNuM 
[SNe (100 yrr^lO^Mo)" 1 ]. 

As with our derivation of the DTD, we will get differ- 
ent rates, depending on the density tracer we use. Table 2 
lists, for each of the density tracers, and for the hybrid com- 
binations, the values or R and Rm in the MCs as a whole. 
Considering the range in R covered by the uncertainties and 
by the different tracers, we see that a SN explodes in the 
MCs once every 200 to 500 years. This is in excellent agree- 
ment with the historical record: the two most recent SNe in 
the MCs were SN1987A (a CC-SN), and SNR J0509.5-6731 
( B0509— 67.5) wh ich exploded circa 1600 and was a SN la 
[ (|Rest et al.ll2008l ); see §5.3 of jBadenes et al.ll2008l )]. It is 
highly implausible that many additional SNe exploded in 
the MCs over the past 400 years but were unnoticed by 
naked-eye observers. 

In terms of mass-normalised SN rates, the range en- 
compassing the different tracers and their la uncertaities, 
1.7 to 3.7 SNuM, is in excellent agreement with rates mea- 
sured in the bluest dwarf galaxies, of which the MCs are 
prototypes. For such galaxi es, which have the hig hest-known 
mass-normalised SN rates, ()Mannucci et al.ll2005l ) have mea- 
sured rates, in SNuM, of 2.3 ± 0.8 (CC-SNe) and 0.9 ± 0.4 
(SNe la), or a total SN rate of 3.2 ± 0.9 SNuM. 

Conversely, we can use the above SN rates to reinforce 
our claim that that SNR sample we compiled in Paper I is 
indeed fairly complete, and that our model, in which the MC 
SNRs are in their Sedov- Taylor phase, is correct. Suppose, 
by way of negation, that the MC SNRs are actually in free 
expansion, as has been invoked by a number of authors to 
explain the observed uniform SNR size distribution (see Pa- 
per I). Since free-expansion SNR shock velocities are in the 
ranee ~ 5000 to 10,000 km s" 1 , the maximum observed SNR 
radii of r ~ 30 pc would translate to maximum ages of 3 to 
6 kyr. The 77 SNRs in the MCs would then imply a SN ex- 
plosion in the MCs once every 40 to 80 years, on average, or 
even more frequently if our SNR sample is incomplete. Sim- 
ilarly, the 3 to 6 kyr visibility time, independent of ambient 
density in this picture, would give mass normalised total SN 
rates of 8 to 15 SNuM, or even higher if our sample is incom- 
plete. This is much higher than in any known type of galaxy. 
Most seriously, perhaps, if the visibility time were fixed at, 
say, 6 kyr, this would raise in the "no scaling" model of 
Table 2 from a value of 2.86 to 6.4 SNe yr" 1 (10 10 M ) _1 , 
corresponding to a CC-SN yield of 0.022 M" 1 . To produce 
so many CC-SNe with any standard IMF, all stars with ini- 
tial mass above 3.7M0 would have to undergo core-collapse. 



This is in direct contradiction to stellar evolution theory, 
and t o the semi-empirical initial-final mass relation for WDs 
(e.g. (I Catalan et al ] |2008l ; ISalari~ et al ] |2009l ; IWilliams" et al.l 

mm) 

In summary, the Sedov model we have adopted to ex- 
plain the MC SNR size distribution is consistent with the 
gas density distributions of the MCs (Paper I), and results 
in MC SN rates that are consistent with the historical record 
and with the rates measured in similar types of galaxies. In 
contrast, a free-expansion model for the MC SNRs produces 
rates that are in stark contradiction with these observations, 
and with expectations from stellar evolution. 



8 SUMMARY AND CONCLUSIONS 

In Paper I, we assembled a multi-wavelength compilation 
of the 77 known SNRs in the MCs, collected from the ex- 
isting literature. We verified that this compilation is fairly 
complete, and that the size distribution of SNRs is approxi- 
mately flat, within the allowed uncertainties, up to a cutoff 
at r ~ 30 pc, as noted by other authors before. We then 
proposed a physical model to explain this size distribution. 
According to our model, most of the SNRs are in their Se- 
dov stage, quickly fading below detection as soon as they 
reach the radiative stage. Under these circumstances, a flat 
distribution of SNR sizes can be obtained, provided that the 
distribution of densities in the MCs follows a power law with 
index — 1. Finally, in Paper I, we used three different density 
tracers (HI column density, Ha flux, and SFR) to demon- 
strate that the distribution of densities in the MCs indeed 
follows a power law of index —1. 

In this paper, we have used the Paper I sample of SNRs 
as an effective supernova survey, conducted over tens of kyr. 
We have applied a novel technique to this SNR sample to 
derive, in these galaxies, the SN rates and the SN DTD. In 
order to accomplish this, we have used the three tracers of 
the density to scale the visibility time of SNRs in the MCs. 
From the scaled visibility times and the SFH maps from 
the resolved stellar pop ulati ons in the MCs publishe d by 
lHarris fc Zaritskvl |2004 ) and lHarris fc Zaritskvl (|2009l ). we 
have calculated the DTD. The DTDs we have derived are the 
first obtained using SFHs from resolved stellar populations. 

Our main findings are the following: 

1. We detect at the > 99% confidence level a "prompt" 
SN la population, defined here as one that explodes within 
330 Myr of star formation. This finding joins a growing num- 
ber of measurements of such a component (see §1). However, 
our measurement is based on our statistically robust DTD 
recovery method, with its avoidance of the averaging inher- 
ent to many previous measurements, and using the most re- 
liable SFHs, based on resolved stellar populations. The yield 
of prompt SN la, in terms of SNe per stellar mass formed, 
is also consistent with other measurements. 

2. We obtain upper limits on the delayed SN la population, 
which are consistent with other measurements. Using these 
upper limits, we find that roughly half, but possibly more, 
of SNe la are "prompt" , as defined above. This again joins 
previous results on the large relative fraction of the prompt 
population. 

3. We use our SNR sample and our SNR visibility times to 
derive current MC SN rates. The SN rate in the MCs agrees 
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well with the historical record for these galaxies. The mass- 
normalised SN rate in the MCs is in excellent agreement 
with the rates measured by SN surveys in galaxies of this 
type. This lends support to the physical model we have pre- 
sented to explain the SNR size distribution and the clouds, 
and which we have used to derive visibility times for our 
DTD and rate calculations. 

4. Conversely, a "free expansion" model for the SNRs, as has 
been invoked for the MCs and other galaxies, would imply 
SN rates in conflict with the historical record and with the 
rates in other star-forming dwarf galaxies. Furthermore, this 
model would indicate an unreasonably high yield of CC-SNe 
per stellar mass formed. 

The main limitation of our study has been the relatively 
small number of SNRs in the MC sample, which has forced 
us to use coarse time bins and has led to large Poisson er- 
rors. Construction of significantly larger samples of SNRs 
is, however, possible by means of deep radio surveys of ad- 
ditional galaxies that are near enough for deriving SFHs 
of their individual regions via resolved stellar populations, 
namely M31 and M33 (see Paper I). Such data would permit 
a similar analysis to the one we have done, but with larger 
SNR numbers, permitting a more accurate determination of 
the SN DTD, and bringing into better focus the properties 
of SNe, their remnants, and the connections between them. 
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